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£N| ■ We have studied the enhancement of the superconducting transition temperature, T c , in a t — 

J — U model of electrons moving on a square lattice in which anisotropic electronic hopping is 
introduced. The inclusion of such hopping mimics, in a approximate fashion, a potentially important 
characteristic of materials possessing stripelike charge and spin correlations. For this model we 
have calculated T c for singlet pairing using the non self-consistent Thouless criterion, and find a 
dramatic enhancement of T c induced by hopping anisotropies. Further, the maximum increase in 
T c is obtained when the system is pushed towards the extreme anisotropy limit, that is, when the 
hopping of electrons is confined to occur in 1 + + dimensions. We demonstrate that in this limit 
the increase in T c , with respect to the isotropic system, can be of the order of 1000. We have also 
determined that in the extreme anisotropy limit the superconducting gap is an equal mixture of s 
^ , and d pairing symmetries (two choices of such a combination being s + d and s + id) owing to the 

S_h ■ reduced (square to rectangular) symmetry of the system in the presence of hopping anisotropies. 

Oh 1 Thus, the presence of d-wave superconducting features in materials whose symmetry is very different 

from that of a two-dimensional square lattice, with the anisotropy produced by the appearance of 
stripes, is not unexpected. 



(N 



C/3 



G 

O 

o 



> 

oo 



s 

I 



X 



I. INTRODUCTION 



Strong support for the existence of stripes in high-T c cuprates, and other transition-metal oxides, has been provided 
by many experiments [jjj. This leads to the question: Do stripe correlations help, hinder, or even possibly create the 
pairing instability that leads to superconductivity? In this paper we demonstrate that at least one feature of stripclikc 
correlations strongly enhances the superconducting transition temperature, T c . 

A previous examination of the magnetic properties of the very weakly doped cuprates || modeled the observed 
experimental support for stripe correlations using an effective Hamiltonian in which a (spatially) anisotropic 
exchange interaction was implemented to represent the stripe-induced magnetic energy scales. That is, in the direction 
parallel to the stripes the full local Cu — Cu exchange would be present, while perpendicular to the stripes a reduced 
exchange would be encountered across such stripes. Renormalized Hamiltonians of a similar simplifying spirit were 
CD also used in other studies of the doped cuprates ■ 

In a recent paper we introduced a model that mimics one aspect of the stripe correlations by incorporating an 
anisotropic hopping Hamiltonian (7); that is, carriers are expected to be able to move much more readily along the 
direction of the stripes, so-called rivers of charge J3|, than perpendicular to the stripes. Allowing the carriers to interact 
via Heisenberg superexchange (using a two-dimensional t — J model), and excluding double occupancy (U — > 00), 
we have found exact solutions for the two electron bound-state problem, viz., in the dilute electron density limit. 
Our calculations demonstrated how stripelike hopping anisotropies can produce (i) a vanishing value for the threshold 
exchange coupling, J c /t, that is required for the stabilization of two-particle bound states, and (ii) a dramatic increase 
in the two-electron binding energies . The physical interpretation of these results is particularly interesting — for 
O ■ the same problem in either one dimension, or in two dimensions with isotropic hopping, a critical superexchange 
of J c /t = 2 is obtained. However, in the extreme anisotropy limit the electrons are confined to hop in only one 
lattice direction, and for such a system we found that an infinitesimal J/t produces bound states. The explanation 
of this result follows from recognizing that in the bound state the two electrons will travel in opposite directions 
on near- neighbour "chains", and when they are on near neighbour sites, and have anti-parallel spins, a S + S~ spin 
exchange interaction flips the chains on which the electrons are travelling. An infinitesimal J/t using this mechanism 
is sufficient to produce a two-electron pair. It is appropriate to think of this result as binding produced by dimensional 
confinement. 

This interesting result naturally leaves us with the task of investigating the transition to a superconducting phase 
in a system with a nonzero electronic density as the system undergoes a crossover from a two-dimensional to a 
l+0 + -dimensional system induced by hopping anisotropies. 

In this paper we present the results of our investigation on the effects of such stripelike hopping anisotropies on 
T c . We use the non self-consistent Thouless criterion S to determine T c , a procedure that is known to reproduce the 
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BCS result for the superconducting transition temperature |B|. We find a dramatic enhancement of T c the greater 
the degree of hopping anisotropy, and the maximum T c obtained is found to saturate in the extreme anisotropy limit. 
For certain system parameters, this enhancement can be of the order of 1000! Thus, this simple model of hopping 
anisotropy, which may be an essential ingredient to stripe correlations, provides a robust demonstration supporting 
the conjecture that stripes can indeed augment pairing correlations. 



II. MODEL HAMILTONIAN 



We model stripelike hopping anisotropy in a square lattice of strongly interacting electrons using the following 
t — J ~ U model: 

H = - %( c L c j> + h - c -) + 53^ij(Ss • Sj - -THTij) + U^2n iA n hl . (1) 

In this Hamiltonian, the sites of a two-dimensional square lattice of size L x x L y with periodic boundary conditions 
are labeled by the indices i and j, tij and Jij are the hopping integrals and exchange couplings between sites i and j, 
respectively, Ci i<T is the annihilation operator for electrons at site i of spin a, rii_ G is the number operator for electrons 
at site i with spin a, and U is the on-site Hubbard energy. 

The most familiar strong-coupling variant of the Hubbard model is the t — J model, and the physics of (square 
lattice) doped Mott insulators described by this model was reviewed by Dagotto (To). As emphasized by, e.g., Anderson 
PI , a vital component of the t — J Hamiltonian is the constraint of no double occupancy. That is, in the t — J model 
one does not use the electron creation and annihilation operators of Eq . ([I]) , but rather one uses constrained creation 
and annihilation operators (for example, see the discussions in Ref. [jTo|). However, the above t — J — U Hamiltonian 
can be used to accomplish this same mathematical projection by taking U — > oo — this simplifying approach has 
been noted by a variety of researchers (see, e.g., Rcfs. and will also be used by us. 

In this paper we restrict ty and to be nonzero for near neighbours (NN) only. Further, as we did in Ref. 0, we 
allow the hopping integral in the x direction, t x , to be different than the hopping integral in the y direction, t y . We 
have also investigated the physics that arises when J x is allowed to be different than J y , but find that no qualitatively 
new physics arises as long as both J x and J y remain nonzero. Thus, from now on we set J = J x — J y , and we analyze 
the resulting hopping anisotropy problem in terms of 

tx = * , r= t -f . (2) 

Thus, we have three dimensionless energy scales in the problem, viz., U/t, J/t, and r; the U/t — » oo limit reduces this 
number to two. 



III. IDENTIFYING T c IN THE LADDER APPROXIMATION 

We use the non self consistent Thouless criterion i§] to identify the temperature at which our system becomes 
unstable with respect to a low temperature superconducting phase. To this end, we determine the equation for the 
two-particle vertex, or effective interaction, in the singlet channel. It is known that at low electron densities |l5| ] 
one may evaluate the vertex function via the ladder approximation to the Bethe-Salpeter equation, and the lack of 
convergence of the sum of the ladder diagrams identifies the critical temperature, viz. when the sum converges the 
system should be in the normal state. The non self consistent formulation of the ladder approximation allows us to 
work exclusively in the normal state, and has been shown to lead to an identical transition temperature to that found 
in the BCS theory of superconductivity ||. 

Thus, we focus on evaluating the effective interaction T when all particle-particle diagrams are included. The 
integral equation for this function in this approximation can be written as (note that we only need to consider r-TT-i-i 
p6| , and to eliminate the proliferation of spin indices in the equations that follow from now on wc suppress all spin 
dependencies from our equations) 

r(fc, k>, Q) = V(k - k') - J -^V(k - q )G°(q)G°(Q - q)T(q, k< , Q) . (3) 

The above integral equation is known as the Bethe-Salpeter equation in the ladder approximation Jl5| , and the various 
functions appearing in this equation are defined as follows: G (q) is the zeroth order single-particle Greens function 
defined by, 
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G°(q) = G°(q,iu; n ) = 



(4) 



the fermionic Matsubara frequencies are given by u) n = (2rt + l)-|; for our anisotropic hopping model the single particle 
dispersion of non-interacting band electrons is given by, 



-2i(cos q x + r cos q y ) ; 



(5) 



/j, is the chemical potential. 

To proceed to the solution of the Bethe-Salpeter equation, we need to reduce it according to symmetries. That is, 
the anisotropy that is present in Eq. ([!]) when < r < 1 reflects a lowering of the point group symmetry of the system 
from that of a square to that of a rectangle. Recall that the basis functions corresponding to the relevant irreducible 
representations of the two-dimensional square lattice are 1, (cos k x + cos k y ) and (cos k x — cos k y ), which correspond to 
on-site s-wave, extended s-wave, and d-wave gap symmetries, respectively. Then note that in our anisotropic model 
these symmetries arc mixed and map onto the fully symmetric A\ irreducible representation of the rectangular point 
group. As a result, it is helpful to decompose the interaction, in the singlet channel (even in k), according to the basis 
functions of A\ of the rectangular point group symmetry. That is, we artificially set J x J y , analyze the resulting 
interaction in terms of a linear combination of the A\ basis functions of 1, cosk X: and cos/cy, and then, in the last 
step of the calculation, restore the square lattice symmetry of the superexchange interaction by resetting J x = J y . 

Following this prescription, and focussing on singlet pairing only, we write the bare interaction term V(k — k') as 



V(k -k') = U- 2J^(k)<Mk') - 2J y ^ y (k)(j) y (k') , 
where <fi x (k) = cosk x and 4> y {k) — cosk y . Substituting Eq. (^|) into Eq. (|j|) we obtain 

r(fc, fc', q) = [u - 2j x <t> x {k)4> x {k l ) - 2j y 4> y (k)4> y (k')] 

<rq iu-2j x Mk)Mq)-2JyM k )Mq)]G Q (q)G Q (Q-qMq,k\Q) 



(6) 



(7) 



(2tt, 

The above equation is an integral equation for the vertex T(q, k 1 ', Q), and demonstrates how the irreducible repre- 
sentations of the square lattice group point are mixed in the rectangular phase. To solve this equation we use a 
standard projection technique (e.g., see Ref. p2|), and setting the total centre-of-mass momentum Q = 0, we obtain 
the following set of equations: 



(1 + ^Xo) - ZJxXx - 2 J y x y 

U~X_x (l ^JxXxx) JyXxy 

UXy — 2J x Xxy (1 — ^JyXyy) 



Co(k) 
C x {k) I = 
C y (k) 



(k) - 2J y x y (t> y (k) 

Xx (1 ^JxXxx^x 

(k)) - 2J y x xy (j) y (k) 

Xy ~ ZJxXxyfix 

(k) (i-2j yXyy <t> y m 



(8) 



We will also discuss results obtained in the no double occupancy limit, U/t — » 00, and in this limit Eq. (ph becomes 



Xo — ZJxXx — ZJyXy 

Xx (1 2>JxXxx*) ^>JyXxy 
Xy ~ %JxXxy (1 — ^JyXyy) 



C (k) 
C*(k) = 

V c„(k) 



Xo -2J x Xx^(k) ~ 2J v Xy(f> y (k) 
Xx (1 ^JxXxx^x (k)) - 2J y x xy <t> y (k) 
Xy ^JxXxy^Px 

(k) [\-2JyXMk)) 



(9) 



The various functions appearing in the above two equations are defined as 

O) = ^£G (q)G (-q)r(q,k) > 



Cx = ^^(g)G (q)G°(-q)r(q,k) 
q 

C « = jf E <M?)G°(q)G°(- q )r( q , k) 



(10) 



3 



The various Q = susceptibilities are functions of /x and T, and are defined by 



Xo 

Xx 
Xy 
Xxx 
Xxy 
Xyy 



l^G°(fc)G°(-fc) , 

k 
k 

i^(k)G°(fc)G (-fc) , (11) 

k 

i^^(k)G°(fc)G°(-fc) , 

k 

I^^(k)G°(fc)G°(-fc) , 

fc 

l^^(k)G (fc)G°(-fc) . 
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To study the pairing instability we determine the temperature at which a divergence of the vertex T(q, fc', Q) occurs. 
Clearly, T(q, fc', Q) depends upon the coefficients Go, C x and C y , defined by Eq. (|l0|), and a singularity in them demands 
the vanishing of the determinant of the coefficient matrix appearing on the left hand side of Eq. (^) ; the vanishing 
of this determinant is thus a simple way in which we can identify the superconducting transition temperature, and 
hence T c is obtained numerically as a function of the chemical potential /i. However, a different and perhaps more 
physical way of displaying our results is to determine T c as a function of the electron density per lattice site, n, and 
in the BCS approximation n at T c is given by 



n( M ,T c ) = 2^ l . . (12) 



e — /3c(ek-£0 -|_ I 

K 

Below we use this latter equation to plot T c vs. n. 

Figure |l| shows our results for U/t = and J/t — 1/3 (similar results are found for other ratios of J/t) with a 
hopping anisotropy of r = 0.1, 0.01, and 0.001. (We have set t = 1 eV as a representative energy; also, this allows 
us to express T c in Kelvin.) As is known from the results corresponding to the isotropic case for U/t = 

the low density (/x ~ —At) region is dominated by (on-site) s-wave pairing, that is to say that the maximum T c is 
obtained for s-wave pairing, whereas near half filling (fi ~ 0), the <i-wave instability dominates. As r deviates from 
the isotropic case, the two pairing symmetries mix and we find a decrease in the values of T c near half filling, and 
an enhancement in T c at low densities. For r = 0.1 we see that there still exists a signature of two broad transitions 
(d-wave and on-site s-wave). As r decreases further the two transitions merge with a T™ ax ~ 200 A'. 

Intermediate ratios of U/t lead to results that smoothly interpolate between U/t — and U/t — oo, and in the 
interests of brevity we only show these two limiting values of U/t. U/t — oo excludes the possibility of an on-site 
s-wave order parameter, hence the pairing symmetry has predominantly extended s- and c?-wave components (see the 
discussion and numerical results in the next section). Our T c vs. n results for these parameters are shown in Fig. ^ 
for the same hopping anisotropies as in Fig. El 

In both of these figures, it should be noted that when we take the hopping anisotropy to be even smaller than 
r = 0.001, no noticeable change occurs — thus, our numerical results for this ratio are representative of the r — > 
extreme hopping anisotropy limit. 

The maximum T c at low densities, viz. T™ ax ~ 70AT (again taking t — leV) occurs at about n — 0.08. As 
a striking demonstration of the effectiveness of stripelike hopping anisotropy in increasing T c , note for isotropic 
hopping, J/t = 1/3, U/t = oo, and a density of n = 0.08 electrons per site, we find that T c ~ 0.08A. That is, we 
find an enhancement of T c due to hopping anisotropy of a factor of about 1000! Clearly, this dramatic increase in 
the superconducting transition temperature supports the conjecture that stripelike correlations can strongly affect 
pairing, and, at the very least, can greatly augment the stability of the superconducting phase at higher temperatures. 
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FIG. 1. The dependence of the superconducting transition temperature, T c , vs. the electronic density per lattice site, n, for 
U/t = and J/t — 1/3. We have set t — 1 eV and express T c in Kelvin. The three curves correspond to anisotropic hopping 
ratios, r, of r = 0.1, 0.01, and 0.001. 




We note that results demonstrating an enhancement of T c due to orthorhombic distortions have been given earlier 
by Li et al. |lS|j , although the extreme anisotropy limit was not considered; the importance of this latter limit was 
made apparent to us in our study of the bound-state formation for the same Hamiltonian as in this paper, albeit for 
only two electrons . 

We also note that the onset of superconducting pairing caused by interchain single particle tunnelling (1ST) has 
been considered earlier by Bourbonnais and Caron pQ| . Starting from a Luttinger Liquid model (a linear array of 
conducting chains separated by some distance) , it was shown in the limit of small interchain hopping 1ST leads to an 
effective pair tunnelling which may eventually induce superconductivity in the singlet channel below a temperature 
T x i (signifying a crossover to a higher dimensionality, viz., a Luttinger Liquid to Fermi Liquid transition). 

In the quasi one-dimensional systems studied in Ref. it was also found the transition temperature, T c , ap- 
proximately scaled as tj_, the interchain hopping integral, (or as (t±) a where a(> 1) is a continuous function of the 
interaction parameter). Our Fig. ^ shows the variation of the maximum T c (as a function of electron density n) vs. 
hopping anisotropy, r, for U/t = 0,1/3 and oo, and for J/t = 1/3. The dependence is found to be roughly linear, 
thus providing some support for the conjecture made by Bourbonnais et al. pCj]. 
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FIG. 3. The maximum T c , found as a function of electronic density for each hopping anisotropy, is plotted as a function of 
the hopping anisotropy r, for U/t = 0, 1/3 and oo and J/t = 1/3. As in Figs. 1 and 2, we have expressed T c in Kelvin using 
t = 1 eV. 



IV. BCS GAP EQUATION 

As mentioned earlier, anisotropic hopping couples the superconducting gaps in the s- and d-wave channels, and 
this leads one to think in terms of a mixing between the two symmetries having the following forms jl9 
s + id state, where there is a phase difference of 7r/2 between the gap functions in s- and d-channels, and (ii) an s + d 
state, where the phase difference is zero. 

We have evaluated the ratio of these gaps by solving the zero temperature BCS gap equation, which can be written 
as, 

A(k) = -£V(k,k')^ , (13) 



where £4. = -y/( £ k A*) 2 + | A.(ki) | 2 . Note that in the U/t — > 00 limit, the on-site s-wave component is suppressed and 
one is left to consider only the extended s- and the <i-wave symmetries. For simplicity, in this section we only report 
on this region of parameter space (U/t — > 00). This allows us to write the trial gap function in k-space with arbitrary 
phase difference 9 in the form, 

A(k)=A s / s (k) + e ie A d / d (k) , (14) 

with A s and A^ representing the amplitudes of the gap functions, which are real, and 

f s (k) = cos k x + cos k y , (15) 
/d(k) = cosfer - cos k y . 

Ignoring the on-site term, we write the interaction in terms of the above-defined basis functions: 

V(k, k') = V(k - k') = -2 J[cos(fc x - k' x ) + cos(k y - k' y )} (16) 
= -J[/ s (k)/ s (k') + / d (k)/ d (k')] • 

Using Eqs. dl^-|l6|) in Eq. (|lj), one can take the real and imaginary parts of Eq. (|l^), and thus obtain coupled 
equations for A s and A^. We have solved these equations numerically for both choices of the relative phase 9, viz. 
9 = ir/2 (corresponding to s + id) and 9 — (corresponding to s + d). 

Our numerical results for the ratio of gap amplitudes for s + id are shown in Fig. ^| — since the results for s + d are 
qualitatively very similar, for brevity we omit that plot. The r = 0.1 case corresponding to the s + d symmetry shows 
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a variation in the ratio A s /Ad from ~ 0.8 (~ 0.6 for s + d) for electronic densities near half filling, to ~ 1.3 (~ 1.3 
for s + d) for fi ~ — 2i(l + r) (low densities). Thus, at low densities the extended s-wave contribution is somewhat 
larger, whereas near half filling the d-wave component is stronger. As the anisotropy is increased, the ratios A s /A<2 
for both s + id and s + d become almost equal to unity for all densities. 
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FIG. 4. The ratio of the gap amplitudes, viz. A 3 /A,j, vs. the chemical potential (in units of t), for an s + id order parameter 
with r — 0.1,0.01 and 0.001. The other parameters are chosen to be U/t = oo and J/t = 1/3. 

We note that similar results have been obtained earlier using variational Monte Carlo studies using a t — J model 
|l9f where the authors found no distinguishable difference between s + d and s + id pairing symmetries. More generally, 
they found that this result was true for any arbitrary 9 in Eq. (|l3|). 

At this stage we cannot ascertain which of the above pairing states is chosen by the system in presence of hop- 
ping anisotropy, so there is a necessity to calculate the condensation energies for these pairing symmetries using a 
microscopically derived Ginzburg-Landau free energy functional — we leave this problem for a future publication. 



V. CONCLUSIONS 



To conclude, we find a dramatic enhancement of the superconducting transition temperature, T c at low electronic 
densities due to stripelike hopping anisotropics. Also, we find that the maximum T c is found to saturate as r — ► 0, that 
is, in the extreme anisotropy limit. In fact, one can show that T c tracks the electronic density of states, consistent with 
the fact that our determination of the superconducting instability is essentially a mean-field theory. As the hopping 
anisotropy mixes different pairing symmetries, such as on-site s-wave, extended s-wave, and d-wave, it is expected 
that the system may choose either a s + d, or s + id, or more generally a pairing of the form A s (k) + e^A^k), 
where 9 may be a function of temperature and anisotropy. To shed some light on this issue we have solved the zero 
temperature BCS gap equation and have evaluated the ratio between the gap amplitudes. Our results demsontrate 
an equal mixing between the s and d gap functions in the extreme anisotropy limit. 

We do not wish to give the impression that we believe that this model is an adequate way of modelling the effects 
that stripelike rivers of charge create. For example, we have ignored all features associated with the discreteness of 
the stripes. Also, although we are including a near-neighbour antiferromagnetic exchange, the ladder approximation 
that we are using to determine a pairing instability will not be adequate near half filling, and thus our theory is not 
just for a moderately doped antiferromagnet. 

Nonetheless, the results that we published in Ref. on two-electron bound state formation, as well as the substantial 
increase in T c that we are presenting in this paper, demonstrate the possible importance of including stripclikc hopping 
anisotropics in realistic microscopic theories. 
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